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Computations  of  ESR  Coupling  Constants  in  Organic  Radicals. 
II.  Nitro  Derivatives  of  the  Phenyl  Radical 


H.  F.  Hameka*  and  F.  J.  Seiler 

Research  Laboratory.  U.S.  Air  Force  Academy.  Colorado  Springs.  CO  80840  (USA) 


Berechnungen  von  ESR  Parameter  in  organischen  Kadikalen.  ll.Nitro- 
phenyl-Radikale 

Wir  beschreihen  die  nach  STO-3G  herechneten  Molekiilstrukturen 
und  ESR-Parameter  fur  zehn  Nitrophenyl-Radikale  und  zwar  die  drei 
Mononitrophenyl-Radikale.  die  sechs  Dinitrophenyl-Radikale  und 
das  2.4.6-1  rinitrophenyl-Radikal  Es  giht  keinc  experimentellen 
Daten. 


Calcul  des  constantes  de  couplage  ESR  dans  des  radicaux  organiques. 
11.  Derives  nitres  du  radical  phenyle 

On  a  determine  par  le  calcul  les  constantes  de  couplage  ESR  el  la 
structure  moleculaire  de  differents  derives  nitres  du  radical  phenyle.  II 
s'agit  en  Poccurence  des  trois  mononitrophenvles.  des  six  dinitrophe- 
nyles  et  du  2,4.6-trinitrophcnvle.  Ces  ealculs  ont  etc  clfeetucs  d'apres 
STO-3G.  II  n  existe  pas  de  resultats  experimentaux. 


Summary 

We  compute  the  ESR  coupling  constants  and  the  geometries  of  a  set 
of  nitro  derivativ  es  of  the  phenyl  radical.  We  consider  the  three  mono- 
nitrophenvls.lhe  six  dinitrophenyls  and  2.4,6-trinitrophenyl.  The  com¬ 
putations  are  based  on  the  use  of  the  Gaussian  82  Program  Package 
with  the  STO-3G  basis  set.  There  are  no  experimental  data  available 


1.  Introduction 

In  the  thermal  decomposition  reactions  of  energetic  mate¬ 
rials  many  of  the  intermediate  products  are  free  radicals,  elec¬ 
tron  spin  resonance  (ESR)  is  therefore  an  effective  and  popu¬ 
lar  technique  for  monitoring  those  reactions.  The  detection  of 
free  radicals  by  means  of  ESR  is  straightforward  hut  the  iden¬ 
tification  of  the  free  radicals  requires  the  knowledge  of  their 
ESR  spectra,  either  from  experiments  or  from  theoretical  pre¬ 
dictions. 

In  the  thermal  decomposition  of  trinitrobenzene  and  trini¬ 
trotoluene  many  of  the  intermediate  products  are  nitro  deriva¬ 
tives  of  the  phenyl  and  benzyl  radicals.  Since  no  experimental 
data  on  the  ESR  spectra  or  on  the  structures  of  these  nitro 
derivatives  are  available  it  would  seem  to  be  useful  to  try  and 
predict  them  from  theoretical  methods. 

We  believe  that  ESR  coupling  constants  may  be  derived 
from  unrestricted  Hartree-Fock  (UHF)  computations  if  the 
latter  are  corrected  hv  removing  spin  contamination.  In  previ¬ 
ous  work11"1  we  presented  a  simple  procedure  for  removing 
the  quartet  spin  contamination  from  a  UHF  wave  function. 
We  applied  the  procedure  to  some  small  organic  radicals  such 
as  methyl,  ethyl,  ally  I  and  vinyl1"  and,  subsequently,  to  the 
phenyl  and  benzyl  radicals  and  to  the  benzene  and  toluene 
radical  anions1". 

Since  the  derivation  and  the  details  of  our  procedure  were 
presented  previously"  "  we  just  describe  the  final  results  here. 
We  make  use  of  the  UHF  method  and  we  assume  that  the 
differences  between  the  a  spin  orbitals  <J> ",  and  the  |3  spin  orbi¬ 
tals  <t>f  are  small  when  i  =  1.2, ...  N  and  that  there  is  one 
unpaired  spin  orbital  We  denote  the  UHF  spin  density 

at  the  position  of  nucleus  Ev  pUHF(a)  and  the  corrected  spin 
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density  by  (yc„r(a).  The  relation  between  t\,„  and  yrHK  is  then 
given  bv 

l\.., (a)  =  (l  '3)‘>i  w(a)  +  (2  3)|<t»s  .  ,(a)l:  (1) 

We  have  computed  the  ESR  coupling  constants  from  Eq.  ( I ) 
and  also  the  molecular  structures  of  a  set  of  nitro  derivatives  of 
the  phenvi  radical  by  using  the  Gaussian  82  Program  Pack¬ 
age1'1.  The  computations  were  based  on  the  UHF  procedure 
and  they  all  made  use  of  the  STO-3G  basis  sets.  The  energies 
and  the  expectation  values  (S;)  of  the  spin  operator  S;  are  all 
listed  in  Table  1 . 

The  derivation  of  Eq.  ( 1 )  is  based  on  the  assumption  that  the 
differences  between  the  a  and  (3  spin  orbitals  are  small,  which 
is  equivalent  with  the  assumption  that  the  spin  contamination 
of  the  UHF  wave  function  is  also  small. 

It  may  be  seen  from  Table  1  that  tins  condition  is  not  satis¬ 
fied  in  the  present  situation  since  the  (S:>  values  are  all  signifi¬ 
cantly  larger  than  the  doublet  value  0.75.  We  feel  nevertheless 
that  our  work  serves  a  useful  purpose  even  though  the  accu¬ 
racy  of  the  results  may  he  open  to  criticism.  First,  we  consider 
not  just  one  system  but  ten  different  nitrophenyl  radicals.  We 
hope  that  our  results  will  make  it  possible  to  detect  the 
similarities  and  differences  between  the  different  radicals.  Sec¬ 
ondly.  our  present  work  represents  the  best  possible  ah  initio 
computations  that  can  be  performed  for  systems  as  large  as  di- 
and  trinitrophenyls.  We  hope  that  our  theoretical  results  are 


Table  1.  Energies  (in  terms  of  hartrees 
Nitrophenyl  Derivatives 

;)  and  < S' >  Values 

of  Various 

E  <S-> 

2-nitrophenvl 

-428.038532 

2.4903 

3-nitrophenyl 

-428.03991 1 

2.4709 

4-nitrophenyl 

-427.981082 

1 .9505 

2.3-dinitrophenyl 

-628.851  125 

3.8029 

2,4-dinitrophenyi 

-628.835  726 

3.4037 

2.5-dinitrophenvl 

-628.836463 

3.3950 

2.6-dinitrophenyl 

-628.834  269 

3.4009 

3.4-dinitrophenyl 

-628.848  297 

3.8051 

3,5-diniirophenvi 

-628.837819 

3.3653 

2.4,ii-uiP''.iopiienvi 
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helpful  in  identifying  the  various  ESR  spectra  since  more  pro  UHF  method.  We  found  thm  our  recults  vccte  more  or  less 

cise  theoretical  results  are  not  available  ut  tliij  liaic  consistent  with  similar  work151,  based  on  the  use  of  the 

Recently  we  also  reported  theoretical  work  on  the  structures  restricted  open-shell  Hartree-Fock  method  (ROHF  method) 
of  the  phenyl  and  benzyl  radicals141,  based  on  the  use  of  the  Our  best  dt.ta  were  derived  from  the  6-31G  basis  set  for  the 
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phenyl  radical  and  from  the  4-31G  basis  set  for  the  benzyl 
radical.  In  the  case  of  the  phenyt  structure  the  6-31G  results 
and  the  STO-3G  results  were  quite  similar.  The  largest  dis¬ 
crepancies  occurred  for  the  CH  bond  lengths,  the  STO-3G 
values  being  about  0.01  angstrom  larger  than  the  6-31G  val¬ 
ues.  However,  the  two  sets  of  bond  angles  were  almost  iden¬ 
tical. 

It  may  be  instructive  also  to  consider  the  theoretical  ener¬ 
gies  that  are  listed  in  Table  1.  In  the  case  of  the  mononi- 
trophenyls  the  meta  derivative  2-nitrophenvl  has  the  lowest 
energy  and  the  para  derivative  4-nitrophenyl  has  the  highest 
energy,  the  difference  is  37  kcal/rnol.  The  energy  of  the  meta 
compound  is  slightly  lower  than  the  energy  of  the  ortho  com¬ 
pound.  the  difference  is  0.9  kcal/rnol.  The  energy  differences 
between  the  various  dinitro  derivatives  are  smaller  than  for  the 
mononitro  derivatives,  the  difference  between  the  lowest 
energy  for  2.3-dinitrophenyl  and  the  highest  energy  for  2,6- 
dinitrophenyl  is  10.6  keal/mo).  Here  the  energies  become 
higher  the  further  the  nitro  groups  are  separated. 

We  present  our  theoretical  results  for  the  ESR  coupling 
constants  in  Figs,  la  to  lj  and  we  discuss  them  in  Sect.  2.  the 
molecular  structure  results  are  shown  in  Figs.  2a  to  2j  and  they 
are  discussed  in  Sect.  3. 

2.  F.SR  Coupling  Constants 

We  computed  the  ESR  coupling  constants  of  the  ten  nitro 
derivatives  of  the  phenyl  radical  by  using  Eq.  ( I ).  The  compu¬ 
tations  are  based  on  the  Gaussian  82  Program  Package' 11  and 
on  the  use  of  the  UHF  procedure  with  the  STO-3G  basis  set.  It 
should  be  noted  that  the  unpaired  orbital  <f>N  .  ,  is  not  neces¬ 
sarily  the  occupied  orbital  with  the  highest  Hartree-Fock 
energy.  We  determined  the  orbital  <f>s  .  i  from  symmetry  con¬ 
siderations  and  from  a  detailed  analysis  of  all  the  a  and  the  (3 
spin  orbitals. 

The  ESR  coupling  constants  were  derived  by  evaluating  the 
values  of  the  unpaired  orbital  .  ,  at  the  positions  of  the 
various  nuclei  by  making  use  of  the  detailed  form  of  the  STO- 
3G  orbitals  and  of  the  definitions  reported  by  Hehre.  Stewart 
and  Pople"’1.  The  coupling  constants  are  then  expressed  in 
terms  of  Gauss  by  using  the  following  conversion  factors: 
1594.04  for  hydrogen.  401  13  for  carbon.  230.43  for  nitrogen 
and  1080.88  for  oxygen.  The  results  are  shown  in  Figs,  la  to  lj. 

We  can  only  speculate  about  the  accuracy  of  each  individual 
value  that  we  report  but  we  hope  to  detect  certain  trends  in  the 
values  and  we  hope  that  these  trends  are  more  meaningful 
than  the  individual  numbers.  The  phenyl  radical  is  formed  by 
removing  one  hvdrogen  atom  from  the  benzene  ring  and  the 
unpaired  electron  remains  localized  around  the  a  carbon  that 
has  lost  its  hydrogen.  Consequently  the  a  carbon  should  have 
a  large  ESR  coupling  constant.  This  is  consistent  with  our 
computations,  which  give  values  between  73  and  85  for  the 
mononitrophenvls  and  values  between  71  and  103  for  the  dini- 
trophenyls.  The  largest  value,  103.  occurs  for  2,3-dini- 
trophenyl. 

We  were  surprised  to  find  that  >•■■*  F.SR  coupling  constants 
of  the  oxygens  in  the  NO?  groups  were  all  very  similar.  With  a 
fpw  exceptions  er..li  pair  of  oxvgens  in  a  NO?  group  has  value*- 
of  around  +27  and  -27  for  each  pair.  Ihe  exceptions  are 
found  in  the  para  NO?  groups  in  the  radicals  with  C;v  sym¬ 
metry,  namely  in  3-mononitrophenyl  and  in  7,4,6-trini- 
trophenyl,  here  the  para  oxygen  coupling  constants  are  all 
around  +5.8.  The  nitrogen  coupling  constants  are  all  small, 
the  largest  values  of  around  5  occur  again  in  the  para  NO? 
groups  of  the  symmetric  radicals. 


It  should  be  noted  that  an  ESR  experiment  does  not  differ¬ 
entiate  between  positive  or  negative  values  of  the  ESR  cou¬ 
pling  constants.  Most  of  the  carbon  coupling  constants  seem  to 
fall  in  the  range  between  15  and  25  if  we  consider  their  abso¬ 
lute  values  only.  Unfortunately,  a  measurement  of  the  carbon 
coupling  constants  does  not  seem  to  be  of  much  help  in  differ¬ 
entiating  between  the  various  nitro  derivatives.  The  excep¬ 
tions  are  the  two  dinitrophenvls  with  adjacent  nitro  groups, 
namely  2,3-  and  3,4-dinitrophenyl.  Here  the  carbon  coupling 
constants  are  a  bit  larger,  varying  between  24  and  33  How¬ 
ever,  if  we  consider  the  accuracy  or  lack  of  accuracy  of  our 
theoretical  results  we  must  admit  that  carbon  coupling  con¬ 
stants  do  not  offer  reliable  guidelines  for  identifying  the  vari¬ 
ous  nitro  derivatives. 

The  coupling  constants  of  the  hydrogens  are  much  easier  to 
measure  than  those  of  the  other  atoms.  The  hydrogen  results 
are  therefore  the  most  important.  Obviously,  an  ESR  experi¬ 
ment  can  easily  determine  the  number  of  hydrogens  in  a  radi¬ 
cal  and  it  serves  to  differentiate  between  a  mono-,  a  di-.  or  a 
trinitrophenyl  even  if  we  do  not  know  the  values  of  the  cou¬ 
pling  constants. 

We  consider  the  mono-,  di-  and  tri-derivatives  separately. 
There  are  only  two  hydrogens  in  2.4.6-trinitropheny!  and  their 
coupling  constants  are  both  equal  to  - 10.  We  did  not  compute 
any  other  trinitro  derivatives. 

In  the  cr.se  of  the  mononitro  derivatives  the  ortho  and  meta 
radicals  have  very  similar  sets  of  coupling  constants,  in  fact 
they  are  so  similar  that  it  becomes  difficult  to  differentiate 
between  the  two  of  them.  The  paia  derivative  has  a  different 
set  of  coupling  constants  and  it  is  also  symmetric  so  that  it 
should  be  easy  to  identify  it. 

In  the  case  of  the  dinitro  derivatives  there  are  sufficient 
differences  between  the  six  sets  of  ESR  spectra  to  allow  iden¬ 
tification.  We  were  unable  to  detect  any  obvious  trends  among 
the  coupling  constants  other  than  the  fact  that  pairs  of  them 
are  equal  in  the  symmetric  systems.  In  three  of  the  radicals 
(2.3;  2.4;  2.61  there  are  two  large  and  one  small  values,  in  one 
radical  (2.5)  there  are  one  large  and  two  small  values,  in  one  of 
them  (3.5)  they  are  all  small  and  in  one  of  them  (3.4)  they  are 
all  large.  The  above  features  show  sufficient  differences  to 
pei  mil  identificaton  of  the  various  dinitro  derivatives. 


3.  Geometries 

In  optimizing  the  geometries  of  the  various  nitrophenyl  radi¬ 
cals  we  imposed  the  condition  that  the  radicals  all  have  a  pla¬ 
nar  structure,  other  than  that  we  varied  all  bond  lengths  and 
bond  angles.  The  results  of  the  geometry  optimizations  are 
presented  in  Figs.  2a  to  2j. 

It  may  be  seen  that  the  ten  geometries  are  fairly  similar,  the 
exception  is  those  radicals  with  C?s  symmetry  and  a  nitro 
group  in  the  paia  position.  In  these  para  nitro  groups  the  C'N 
bond  length  is  1.506  angstrom  and  the  NO  bond  lengths  are 
1.280  angstrom.  In  all  other  nitro  groups  the  CN  bond  lengths 
are  between  1.46  and  1.47  angstrom  and  the  two  NO  bond 
lengths  differ  by  about  0.015  angstrom,  the  larger  one  being 
phnut  t  V,  -?r>d  the  smaller  one  being  about  1.345  angstrom 
The  CH  bund  lengths  me  leinaikUjl)  constant  in  ail  radicals, 
they  vary  between  1.081  and  1.083  angstrom. 

The  bond  angles  of  the  phenyl  ring  are  not  strongly  affected 
by  the  nitro  groups.  The  a  carbon  usually  has  the  largest  bond 
angle,  the  average  value  for  the  ten  radicals  is  123.7  degrees. 
The  bond  angle  at  the  para  position  is  usually  a  bit  larger  than 
120  degrees.  The  largest  differences  occur  for  trinitrophenyl 
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where  the  a  carbon  has  a  bond  angle  of  122. U  and  the  para 
carbon  has  a  bond  angle  of  123.6  degrees.  The  ortho  and  meta 
bond  angles  are  generally  a  little  smaller  than  120  degrees.  The 
presence  of  a  nitro  group  seems  to  increase  the  bond  angle  bv 
an  amount  between  0.5  and  1.0  degree  as  compared  to  the 
average  value. 


4.  Final  Remarks 

We  can  only  speculate  about  the  accuracy  of  our  computa¬ 
tional  results  but  we  hope  that  they  may  be  helpful  in  detecting 
and  identifying  the  presence  of  nitrophenyl  radical  in  organic 
reactions  and  that  they  describe  at  least  the  general  trends  of 
the  ESR  spectra  and  the  geometries  of  those  radicals. 
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